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Universality in Transport Processes of Unconventional Superconductors 
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We show that some of the low temperature transport coefficients (e.g., electrical and thermal 
conductivities, viscosity and sound attenuation) are universal, i.e., independent of the impurity con- 
centration and phase shift for specific classes of unconventional superconductors. The existence of a 
universal limit depends on the symmetry of the order parameter and is achieved at low temperatures 
feT < 7 < Ao, where 7 is the bandwidth of the impurity induced Andreev bound states. The 
density of states is finite at zero energy and leads to the re-appearance of the Wiedemann-Franz law 
deep in the superconducting phase for fcsT <C 7. Our findings also show that impurity concentration 
studies at low temperatures can distinguish between different order parameter symmetries. 



PACS: 74.25.Fy, 74.25.Ld, 74.72.-b., 74.70.Tx 

We investigate the behavior of the heat current, the 
electrical current, and the momentum current (transverse 
sound attenuation) for an unconventional superconduc- 
tor, i.e., for an order parameter with reduced symmetry 
for which gapless excitations exist even at zero temper- 
ature. Such superconducting states have been argued to 
both exist in the cuprates and heavy fermion supercon- 
ductors. A leading candidate in the cuprates for a tetrag- 
onal crystal structure (D 4 ?j) is the Bi g state (d x 2- y 2), a 
singlet state with lines of nodes at the Fermi positions 
Pf x = ip/y ■ Similarly, the most promising candidates 
in the heavy fermion metal UPt3 are the two-dimensional 
orbital representations coupled to a symmetry breaking 
field. For UPt3, which has a hexagonal crystal structure 
{D 6 h), phase diagram studies, and transport measure- 
ments lead to cither an even-parity, spin singlet Ei g , or 
an odd-parity, spin triplet Ei u pairing state. In both 
cases the order parameter vanishes at the Fermi surface 
on a line in the basal plane, pf z = 0, and at points at the 
poles, pf x =p fy = Q g. 

Instead of examining the effects of the multi-sheeted 
Fermi surface on the heat current, charge current, and 
momentum current, we model the excitation spectrum 
by an excitation gap that opens at line and point nodes 
on the Fermi surface, and by the Fermi surface properties 
in the vicinity of the nodes (i.e., the Fermi velocities, v/, 
and the density of states, Nf, near the nodes) ||. Crystal 
symmetry determines the positions of the nodal regions 
of the excitation gap on the Fermi surface, but not the 
prcfactors for the gap opening. 

The low-temperature behavior of the transport coef- 
ficients probes lower- dimensional regions of the Fermi 
surface, e < 7 <C Aq <C Ef, where the excitation gap 



vanishes, and is less sensitive to the overall geometry of 
the Fermi surface. The low-energy scale 7 is defined by 
the bandwidth of the impurity induced Andreev bound 
states, and reflects the formation of a novel metallic state 
deep in the superconducting phase; for strong scattering 
7 cx yToAo, where Tq is the (elastic) normal-state scat- 
tering rate H/2t(0). We parametrize the nodal regions of 
the gap with a minimal set of nodal parameters, and at- 
tempt to fit these parameters in order to achieve accurate 
low temperature limits for the different transport coeffi- 
cients along the principal axes of the crystal. Thus, our 
order parameter model depends on angle (p/) and the 
nodal parameters ({/Xi}), A(p^;{/ij}). The advantage of 
this approach is that we can quantitatively determine the 
phase space contributing to the low temperature trans- 
port coefficients and then examine in more detail the ef- 
fects of impurity scattering and order parameter symme- 
try on the current response ||, without having to know 
the overall shape of the Fermi surface or basis functions. 

The number of nodal parameters is fixed by the mini- 
mal number of symmetry unrelated point and line nodes. 
These parameters define the slope or curvature of the gap 
near a line or point node in a spherical coordinate system 
(uniaxial anisotropy is included by mapping an ellipsoidal 
Fermi surface onto a sphere). Specifically, for the E-rep 
models of UPt3 we parametrize the gap in the vicinity of 
the equatorial line node by |A(0)| « ^Ao|? — 0\, while 
near the poles |A(0)| ~ /i n Ao|#|™, where the internal 
phase winding number is n — 1 for E\ g and n = 2 for E 2u ■ 
This parametrization allows us to adjust independently 
the opening of the gap at the line and point nodes in or- 
der to describe the low-energy excitation spectrum. Note 
that the crucial difference between the E\ g and Em state 
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lies in the opening of the gap at the polar point nodes. A 
similar parametrization yields for the d x 2_ y 2-w&ve order 
parameter model in the cuprates with a cylindrical Fermi 
surface: A (if) w ^A (j — f), if « ir/A. 

We calculated in linear response (including vertex cor- 
rections), and in the long wavelength limit (q — > 0), 
for sufficiently low temperatures and external frequen- 
cies (fcsT, hui <C 7 <C Ao) the transport coefficients in 
the limit of weak and strong impurity scattering. The 
asymptotic values, which we derived for the electrical and 
thermal conductivity, as well as for the transverse viscos- 
ity rj and ultrasound attenuation a, are listed in Table I. 
Here, we restricted ourselves to the hydrodynamic limit 
and orientations of the wavevector q and polarization e, 
such that atij = (q 2 /gc s )r]ij t ijiiqj 1 where g is the mass 
density and c s the speed of sound. 

Table I. Asymptotic low-T and low-o; values of the transport 
coefficients for three different pairing states, where we have 
used ctoo = e 2 NfVfToo, koo = (ir 2 /3)k%T Nfvfroo, and a o = 
(q 2 1 'Ag c s )Pf NfV 2 Too with an effective too = ?i/2/iAo(0). 
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* In the strong scattering limit including vertex corrections. 

The transport coefficients in the basal plane (or CuC>2 
planes) are universal, i.e., independent of the impurity 
concentration and scattering phase shift. Furthermore, 
the in-plane results do not distinguish between the two 
E-rep models in UPt3. Estimates of these coefficients 
are in good agreement with experiments ||. However, 
transport measurements along the crystal c axis can dis- 
tinguish between an E\ g or Em pairing state. For the 
electrical and thermal conductivity the Em state leads to 
a universal value at low temperatures, independent of im- 
purity scattering, while the E\ g state has a non-universal 
value, strongly dependent on impurity scattering. 

Independent of the universality of the transport co- 
efficients we find that the Wiedemann-Franz law is re- 
established at low temperatures fc^T < 7, due to the 



finite density of states (DOS) at zero energy. At finite 
temperatures the Lorenz ratio L(T) deviates significantly 
from Sommerfeld's value L$ — ^-(fcg/e) 2 even for elastic 
scattering, because of the different coherence factors in 
the electrical and thermal conductivity, and the energy 
dependence of the DOS. The Lorenz ratio depends sensi- 
tively on the scattering phase shift. For strong scattering 
the ratio L(T)/Ls is typically larger than one, whereas 
for weak scattering (Born) it is less than one || . This is 
a very robust feature, and remains true when we include 
inelastic scattering, modeled by a phcnomcnological tem- 
perature dependent relaxation time, as shown in Fig. 1. 
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Fig. 1 Lorenz ratio L(T) for an E2 U state for weak (Born) 
and strong (unitarity) scattering and various phenomeno- 
logical scattering rates L(T) = r (l + T 2 /T 2 ). 

As a result, measurements of transport coefficients 
along different crystal axes at very low temperatures and 
for various impurity concentrations will distinguish be- 
tween the E\ g and Em pairing models in UPt3, and also 
elucidate the pairing state in the cuprates. 



ACKNOWLEDGMENTS 

This research was supported by the NSF through the 
Science and Technology Center for Superconductivity 
(DMR 91-20000), and the Northwestern University Ma- 
terials Science Center (DMR 91-20521). 



[1] D. J. Scalapino, Physics Report 250 (1995) 329. 
[2] J. A. Sauls, Advances in Physics 43 (1994) 113. 
[3] M. J. Graf el at, J. of Low Temp. Phys. 102 (1996) 367; 
Phys. Rev. B 53, 15147 (1996). 



1006 



